arXiv:1507.08596vl [math.DS] 30Jul2015 


Sliding Hopf bifurcation in interval systems 

E. Hooton^’* *, Z. Balanov^, W. Krawcewicz^, D. Rachinskii^ 


D 

D 


Abstract 

In this paper, the equivariant degree theory is used to analyze the occurrence of the Hopf 
bifurcation under effectively verifiable mild conditions. We combine the abstract result with 
standard interval polynomial techniques based on Kharitonov’s theorem to show the existence 
of a branch of periodic solutions emanating from the equilibrium in the settings relevant to 
robust control. The results are illustrated with a number of examples. 


1 Introduction 


Subject and goal. Many problems in population dynamics, neural networks, fluid dynamics, 
solid mechanics, elasticity, chemistry, mechanical and electrical engineering lead to studying the so- 
called Hopf bifurcation (more precisely, Poincare-Andronov-Hopf bifurcation) in dynamical systems 
parameterized by a real parameter (see, for example, [^[^[^1^1^ and references therein). To 
be more specihc, given a parameterized family 


X = f{a,x), a G [a_,a+], x G 


( 1 ) 


where / : [a_, a+] x —>■ is a continuous map and (a, 0) is a curve of trivial stationary solutions, 
the Hopf bifurcation is a phenomenon occurring when a crosses some critical value Uo (for which 
the linearization D^fia, 0) admits a purely imaginary eigenvalue) and resulting in appearance of a 
branch of small amplitude periodic solutions near the curve (a, 0). In his original work [23] , E. Hopf 
studied system Q under the following assumptions: (a) / is analytic in both variables; (b) for 
a = tto, exactly two complex conjugate characteristic roots and intersect the imaginary 
axis (absence of multiple/resonant roots); (c) ^ 0 (exclusion of steady-state bifurcation); and, 

(d) Rey(O) ^ 0 (transversality). Hopf’s theorem includes conditions for the occurrence of the 
bifurcation (i.e., the existence result) and conditions for stability of small cycles bifurcating from 
the stationary point. After this pioneering work, a substantial effort was made in order to relax 
conditions (a)-(d) (see, for example, [^ [I^[^[^|34[|36[|39| and references therein). One objective 
of this paper is to present an abstract result on the occurrence of the Hopf bifurcation in Q under 
very mild (and effectively verifiable) hypotheses containing many known occurrence results as a 
particular case (cf. Theorem |3.2[ Theorem 3.5 and Remark 3.6 1 . It should be stressed that we do 
not study stability of bifurcating periodic solutions. 

Our choice of the conditions on the nonlinearity / and its derivative D^fia, 0) is essentially 
determined by the following observations. In analysis and design, it is customary to deal with 
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approximations of complex models that have some degree of uncertainty (one can think of the so- 
called nominal systems widely used in robust control; see, for example, [^). Considering a model 
with uncertain parameters, one can expect that the entries of the matrix Dxf{a, 0) belong to some 
known intervals of values rather than being represented by fixed numbers. This suggests to study 
the Hopf bifurcation phenomenon for a class of systems Q where coefficients of the linearization 
are limited to known intervals. In this setting, the characteristic polynomial of Dxf{a, 0) that 
defines the stability properties of the linearization also becomes an interval polynomial (see, for 
example, (zl)- Importantly, this setting includes the scenario when the characteristic values of the 
linearization of a representative system Q slide along the imaginary axis when the bifurcation 
parameter is varied (see Figure [^). The main goal of the present paper is to propose a method 
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Figure 1: (a) Given an a-parametrized family of characteristic polynomials with unknown coef¬ 
ficients that are limited to some intervals, the dashed lines bound a corridor for the real part 
T{a) = Re^(a) of an eigenvalue, while the solid line indicates a sliding scenario for some selector 
of the family, (b) Possible complex behavior of the branch of periodic solutions. 


for analysis of the occurrence of the Hopf bifurcation in the presence of such sliding. 

As a matter of fact, the sliding phenomenon makes the problem non-local. Namely, it does not 
allow one to localize a bifurcation point on the basis of the knowledge of the linearization, that is 
based on the condition Re^(a) = 0 (see Figure a,b). 

To study the Hopf bifurcation in this setting, one needs to deal with the whole interval of sliding 
that consists of potential bifurcation points. Thus, sliding is in sharp contrast to the transversality 
condition (d) above. At the same time, to the best of our knowledge, all the existing results on 
the occurrence of Hopf bifurcation identify explicitly a critical value of the parameter a at which 
Re/r(Q;) changes its sign (for the least restrictive condition of this type, we refer to [^). Some 
conditions for the existence of a branch of cycles that are non-local with respect to the parameter 
can be found in [29||^. 

The simplest scenario which includes sliding and is covered by our results is the following. 
Suppose that system 0 has an equilibrium a: = 0 for all values of the parameter a G [a_,a_|_]. 
Assume that the linearization Dxf{ct, 0) of the right hand side is invertible and has at most one pair 
of purely imaginary eigenvalues for any a € [a-,a+\. Finally, assume that the zero equilibrium 
is hyperbolic for a = a± and the dimension of the stable manifold of the linearization of Q at 
zero is different for a = a_ and a = a+. Then there is a Hopf bifurcation point on the interval 
(a_, a+). Theorem 3.2 presented below also covers more complex scenarios including multiple and 
resonant eigenvalues of the linearization on the imaginary axis. 

23 , the Hopf bifurcation in Q was studied based on the series expansion of /. 


Method. In 


The further progress was related to the methods rooted in the singularity theory: assuming that 
the system satisfies several regularity and genericity conditions, one can combine the normal form 
classification with Center Manifold Theorem/averaging method/Lyapunov-Schmidt reduction. For 
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a detailed exposition of these concepts and related techniques, we refer to 20 21 39 


Being very effective in the settings they are usually applied to, the singularity theory based 
methods meet difficulties if a setting is not regular/generic enough. For example, dynamical 
systems with hysteresis components admit linearization at the origin while any small neighborhood 
of the origin contains non-differentiability points which makes the Center Manifold Reduction 
impossible (see [^|^|^[^[^[^|^ ) for details). As long as the stability of bifurcating solutions 
is not questioned, one can use homotopy theory based methods. Important steps in this direction 
were done in (framed bordism theory), [13| (Fuller index), [36| (parameter functionalization 
method combined with the Leray-Schauder degree), to mention a few. 

During the last twenty years the equivariant degree theory emerged in non-linear analysis (for the 
detailed exposition of this theory, including historical remarks, we refer to recent monographs 026 


and surveys 25 


for the prototypal invariants, see 14 15,18 ^). The equivariant degree, 
being the main topological tool used in this paper, is an instrument that allows “counting” orbits of 
solutions to symmetric equations in the same way as the usual Brouwer degree does, but according 
to their symmetry properties. In particular, the equivariant degree theory has all the attributes 
allowing its application in non-smooth and non-generic equivariant settings related to equivariant 
dynamical systems having, in general, infinite dimensional phase spaces with lack of linear structure 
(cf. [^). We refer to and references therein for the equivariant degree treatment of the 

(symmetric) Hopf bifurcation in different environments (see also [28| ). In the present paper, we 
use the A^-degree with one free parameter (see for the axiomatic approach). 

Theorem |3.7| below explicitly refers to the verification of stability properties of interval poly¬ 
nomials (cf. conditions (R3) and (R4)). Among very few results on the connection between 
perturbations of the coefficient and root locations, Kharitonov’s theorem ( [^, see also [7|[^) 
takes a firm position. To be more specific, V. L. Kharitonov showed that given a family of interval 
polynomials with real coefficients, it is necessary and sufficient to test just four canonically defined 
members of the family in order to decide that all polynomials are Hurwitz stable. The main topo¬ 
logical ingredient of Kharitonov’s proof is the so-called Zero Exclusion Principle (in short ZEP) 
which can be traced back to the classical Argument principle in Complex Analysis. In this paper, 
combining ZEP with simple combinatorial arguments, we establish a Kharitonov type result for the 
so-called fc-stable interval polynomials (cf. Lemma 2.4 and Definition 2.2). In particular, it shows 
that Kharitonov’s approach is sensitive not only to Hurwitz stability, but also to the change of the 
dimension of the stable manifold in families of interval polynomials which is crucial for studying 
the Hopf bifurcation phenomenon. 

The paper is organized as follows. In the next section, we present some background related 
to the Hopf bifurcation and interval polynomials. In Section main results are formulated (see 


Theorems 3.2 3.5 and 3. 7|. Some examples illustrating Theorems 3.5 and 3.7 are given in Section 


1^ Section contains the proof of Theorem |3.2| which is close in spirit to the proofs of Theorems 
9.18 and 9.24 from [^. In Section § we provide the proofs of remaining results. A brief summary 
of properties of the (^-equivariant degree is presented in Appendix. 


2 Preliminaries 


2.1 Hopf bifurcation 


The Hopf bfurcation being the main subject of the present paper is formalized in the following 
definition (cf. [5 36 ). 


Definition 2.1. Consider a non-empty set F of non-constant periodic solutions {a.,p,x(f)) of 
system 0 (where p is the minimal period of x{t)) such that p G [p_,p+] C (0, oo). The set F is 
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called a branch bifurcating from the trivial solution if the union ofT and the set of trivial solutions, 
rlj [q;_,q;_|_] X [p_,p_|_] X {x = 0}, is a connected compact set. 


If r is a branch of non-constant periodic solutions bifurcating from the trivial solution, then 
the interval \a-,a+\ contains at least one Hopf bifurcation point ao in the weak sense of [36] . In 
other words, there are converging sequences ak —t ao and pk ^ Po > 0 such that system (1^ with 
a = ak has a non-constant periodic solution Xk{t) with the minimal period pk and Hxfcljc —t 0. 
If the necessary condition for the Hopf bifurcation (see, Section 5.11 is satisfied at exactly one 
point ao G then Definition |2.1| reduces to the definition of the Hopf bifurcation used 

in p. 260]. However, the setting of Definition 2.1 does not exclude a possibility of more complex 
behavior of the branch shown in Figure in the case of an eigenvalue sliding along the imaginary 
axis as in Figure]^. 


2.2 Interval polynomials and Kharitonov’s theorem 

Following [^, an interval matrix, denoted 

= {Ikj)k,j—0 ~ ■ -^kj € Ikj^ h, j = 1 ,..., n} , 

is the set of all matrices whose (fc,j)-th entry lies in the interval Ikj- Similarly, for interval 
polynomials, 

S = Iq IiX ^ InX^ = {H = ao-|-aiA-|-’’’-|- anX^ a^ & Ik, k = , n}. 

Naturally, -1- ^2 = {Pi -I- P 2 : Pi G 5i, P 2 G 52 } and iSi • 1 S 2 = {P 1 P 2 : Pi G iSi, P 2 G 1 S 2 }. 
Also, we will need the following definition. 

Definition 2.2. (i) A polynomial with real coefficients is called monic if its highest coefficient is 
one. An interval polynomial S is called monic if any P € S is monic. 

(ii) Let S be a monic interval polynomial of degree n. We say that S is q-stable (resp., q- 
unstable) if for any P G S, P has exactly q roots with Re (z) < 0 and n — q roots with Re {z) > 0 
(resp., q roots with Re (z) > 0 and n — q roots with Re (z) < 0). 

The classical Hurwitz stability is, therefore, called 0-instability in our terminology. Given an 
interval polynomial 

S = Iq + IiX -I- • • • -I- /„_iA" ^ -|- A", Ij = [aj, bj], (2) 


we denote 


gi{S, A) — cto T b2X? <14A^ 
hi[S, A) = oiA -f 63A^ -I- dsA*^ -!-•••; 
Notice that for any P G S, 


92 ( 5 , A) — &o + U2A^ -I- 64 A"^ -!-•••; 
h2{S, A) = &iA -f a^X^ b^X^ 


Re {gi{S, iuj)) < Re (P(fw)) < Re {g 2 {S, ioj)), Im {hi{S, iur)) < Im {P{iu})) < Im {h 2 {S, iiv)). 

( 4 ) 

The following classical result regarding stability of interval polynomials is known as Kharitonov’s 
theorem (see [7 

Theorem 2.3 (Kharitonov). The interval polynomial ([^ is Hurwitz stable if and only if the 
following polynomials are Hurwitz stable: 


9i(.^,') + '), 9i{S, ■) + h2{S, ■), g2iS, ■) + hi{S, ■), g2{S,-) + h2{S,-). 
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We will use a g-unstable variant of Kharitonov’s theorem. 

Lemma 2.4. If a polynomial Po G S is q-unstable and 

{Re{gi{S,iuj)) ■ Re ( 52 ( 5 , iw)) , Im iw)) • Im (h 2 ( 5 , iw))) ^ {{x,y) ■■ x <0,y <0} 

for any lu > 0, then the interval polynomial S is q-unstable. 

The main topological ingredient of the proof of both statements is the so-called 

Zero Exclusion Principle. If some polynomial Pq G S is g-unstable and for any P G S and any 
w > 0, P(iuj) 7 ^ 0, then the interval polynomial S is g-unstable. 

Proof of Lemma |2.4| As an immediate consequence of inequalities Q one has that, for any 
PgS, 

Re{gi{S,iuj)) ■Re{g 2 {S,iuj)) > 0 ReP{iu;)f^0, 

Im iw)) • Im(h 2 (‘ 5 , lo;)) > 0 => ImP(jtLi) 7 ^ 0 . 

The result then follows from the Zero Exclusion Principle. □ 

2.3 Interval polynomials and Descartes’ Criterion 

Recall the following classical result. 

Descartes’ criterion. If the terms of a single-variable polynomial with real coefficients are ordered 
by descending variable exponent, then the number of positive roots of the polynomial is less than 
or equal to the number of sign differences between consecutive nonzero coefficients. 

As an immediate consequence, we have 

Proposition 2.5. Given a polynomial P with real coefficients, assume that there exist polynomials 
Q and R such that the coefficients of the polynomial 

S{P, Q, R){uj) = Q{uj)Re{P{iu!)) -\- R{u})lin{P(iuj)) (6) 

have at most one sign change. Then, P may have at most one pair of purely imaginary roots. 

Indeed, for w > 0, if iuj is a root of P, then w is a (positive) root of S{P, Q, R). 

In what follows, we use an interval polynomial variant of Proposition |2.5[ For the precise 
formulation, we need the following definition. Given an interval polynomial S, we say that the 
coefficients of S have at most one sign change if, for some j, either [a*,, 6 ^] C (— 00 ,0] for all k < j 
and [ak,bk] C [0, 00 ) for all k > j, or [ak,bk] C (—oo,0] for all k > j and [ak,bk] C [0, 00 ) for all 
k < j. Notice that if the coefficients of S have at most one sign change then the coefficients of any 
polynomial P G S have at most one sign change. 

Set 

T{S, Q, i?)(u;) = Q(a;)Re(5(ia;)) -|- i?(w)Im(5(*a;)). 

Lemma 2.6. Assume that there exist Q,R such that the coefficients ofT{S,Q,R) have at most 
one sign change. Then, any polynomial P G S has at most one pair of purely imaginary roots. 

Proof: Suppose, for the contrary, that some P G S has more than one pair of purely imaginary 

roots. By Q, S{P,Q, R){uj) G T{S,Q,R). Therefore, S[P,Q,R){uj) has at least two distinct 
positive real roots. Hence, by Descartes’ criterion, the coefficients of S{P, Q, R){u!) have more than 
one sign change, which is a contradiction. □ 
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3 Main results 


3.1 Abstract result 

Set ]/ = and assume that / : [a_,a_|_] x F —>■ F is a map satisfying the following properties: 
(PO) / is continuous; 

(PI) The Jacobi matrix D^f^a, 0) exists for all a, depends continuously on a and 


hm sup -- - = 0 ; 

lkl |-!-0 a \x\ 


(7) 


(P2) /(a,0) = 0 for all a; 

(P3) det(Z?a;/(a, 0)) ^ 0 for all a. 

To formulate the next condition, take the map A : [Q!_,a+] x M x K —>■ C determined by the 
characteristic polynomial of the Jacobi matrix Dj;f(a,0), i.e. 

A(a, P, t) = detc(('r + */3)Id - 0)). (8) 

Define a Z 2 -action on by 

(a,/3,r) !-)> (q!,-/3,t). 

Also, given a set P C [«_,«+] x K x 1R+, define 

P± = P Pi ({q;±} X M X IR+) and Pq = P P ([q:_, a+] x M x {0}), (9) 

where K+ denotes the non-negative semi-axis. We will denote by i9fl the boundary of a domain 12 
and by 17 the closure of 17. 

(P4) There exists a bounded Z 2 -invariant domain P C [a_,Q!+] x K x K+ such that: 

(i) P is homeomorphic to a closed ball; 

(ii) A(a, I3,t)^Q for all (a, /3, t) &dV\ (P+ Po); 

(hi) P+ and P_ contain a different number of roots of A(q:, /3, r) (counted according to their 
multiplicities). 

(P5) There exists a finite collection of disjoint sets C [«_, a+] x M_|_ such that: 

(i) each P*. is homeomorphic to a closed disk; 

(ii) A-i(0)n(U77fe X {0}) = A-i(0)nPo; 

(hi) for any I G N and for any {a, (3) G dVk, A(q!,//3,0) ^ 0. 


Remark 3.1. Conditions (PO) and (PI) reflect the minimal regularity that we require from 
system 0 . Condition (P2) guarantees the existence of a branch of zero equilibria from which we 
expect the occurrence of the Hopf bifurcation, while (P3) excludes steady-state bifurcation. 

The domain P provided by (P4) acts as a “trap” catching the roots of A, which may potentially 
contribute to the Hopf bifurcation. Condition (P4)(ii) guarantees that the roots may only escape 
P through the planes {a = a_}, {a = a+} and {r = 0}. Condition (P4)(iii) is an analog of the 
standard non-zero crossing number assumption. 

On the other hand, the sets P/j provided by (P5) form the domain on which we will compute 
the topological invariant. Property (P5)(iii) (which is a kind of non-resonance condition) ensures 
that the topological invariant is well-defined, while (P5)(ii) (which says that all the roots in P^ 
are precisely those “exiting” P) ensures that the invariant is non-trivial and thus that the Hopf 
bifurcation takes place. Several versions of conditions (P4) and (P5) directly related to the 
classical setting for the Hopf bifurcation are discussed in the next subsection. 
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The following statement is our main abstract result. 

Theorem 3.2. Let f satisfy conditions (PO) - (P5). Then, there exists a branch of non-constant 
periodic solutions to system 0 bifurcating from the trivial solution (cf. Definition\2. 


3.2 Corollaries 


Let us consider some corollaries of Theorem 3.2 based on variations of conditions (P4) and (P5) 
which are more relaxed but easier to verify. To this end, we introduce the following notation: 


R{f) = {(a>/3) S [a!_,a+] x M : A(q!,^,0) = 0}, 

Sjif) = {{a, 13) e [a_,a+] x K : (a, j/3) G R{f )}, j = 2,3,..., 

oo (10) 

S{f) = U sAf)- 

t=2 


Remark 3.3. Notice that R{f) is the set of purely imaginary characteristic roots lying between 
a_ and a+, while S{f) is the set of points an integer multiple of which lies in R{f). 

We use a few variants of conditions (P4) and (P5). 

(P4') There exist for which x = 0 is a hyperbolic equilibrium of 0 and the dimension of 

the unstable manifold of the linearization of 0 at 0 is different for a- and a+. 

(P5') There exists a finite collection of disjoint sets Vk C [a_,a+] x IR_|_ such that: 

(i) each Vk is homeomorphic to a closed disk; 

(ii) R{f) c[jPk-, 

(hi) S{f) fl dVk = 0 for any k. 

(P5") R{f) n S{f) = 0. 

(P5'") Dxf{ct,0) has at most one pair of purely imaginary eigenvalues for all a G [Q:_,a+]. 
(P5"") There exists a unique a G (a_, a+) such that Dj,f[a, 0) has purely imaginary eigenvalues. 


Remark 3.4. Observe that (P4') is a non-zero crossing number condition; in particular, the 
classical Routh-Hurwitz criterion (see, for example, |^) can be useful for its verihcation. Condition 
(P5') is a slight modification of (P5), adjusted to the case when (P4') holds. Condition (P5") 
is the classical non-resonance condition. Condition (P5"'), although much more restrictive than 
condition (P5"), can be verified using Descartes’ criterion (see also Proposition |2.5[ ). Finally, 
(P5"") is the standard isolated center condition (see, for example, E)- 

The following statement is based on Theorem 3.2 and is used below to obtain sufficient condi¬ 
tions for the Hopf bifurcation in interval systems. 


Theorem 3.5. Suppose f satisfies conditions (PO) - (P3). Suppose, in addition, f satisfies one 
of the following assumptions: 

(a) (P4') and (P5'); 

(h) (P4') and (P5"); 

(c) (P4') and (P5"'); 

(d) (P4) and (P5"); 

(e) (P4) and (P5"")- 

Then, system 0 has a branch of non-constant periodic solutions bifurcating from the trivial one. 
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Remark 3.6. Under the assumption that / is of class C^’^, Theorem 3.5 ’e) was established in 
(see also mmy On the other hand, by taking a sufficiently small neighborhood (a_, a+), one 
can deduce the main result of [36| from Theorem |3.5K d) (without extra “simplicity” assumptions 
on the corresponding eigenvalues). 


3.3 Theorem 3.5 and interval polynomials 


In this section, we address families of one-parameter systems for which every member is undergoing 
the Hopf bifurcation. To be more precise, denote hy £/ a map from [a_, a+] to the set of interval 
matrices of size d x d and by r a set of maps r : [a_, a+] x V —>■ V. By the symbol 


X = £/(a)x + x) (11) 

we mean the family of all systems of the form 

X = A{a)x + r{a,x) (12) 


satisfying the following conditions: 

(i) A : [a_,Q!+] — L{d,R) is continuous; 

(ii) A{a) G sA (a) for every a G M; 

(iii) r G V. 

Denote by O the map from M to the set of monic interval polynomials such that for any a G K, 
£l{a) = Jo{a) + Jiia)x-\ -h J„_i(a)a;"“^-|-x”, Jfc(a) = [afc(a), 6fc(a)], (13) 


is the collection of all possible characteristic polynomials corresponding to each member of the 
family £/{a) (in fact, this collection constitutes an interval polynomial). To generalize Theorem 

Given a family of 


3.5 a,b,c) to the interval setting, we need “interval analogs” of notations ([10| 
systems ([TI|) with interval characteristic equation (13), put (cf. ^ and ([^) 


wi{a,l3) = Re(5i(0(a),i/3)) • Re(g2(0(a),i/3)), 

W2{a,l3) = Im(/ii(0(a),i/3)) • Im(/i2(0(a),i/3)), 
m{£/) = {(a,/3) : a G [a_, a+], uii(a,/3) < 0,W2{a,l3) < 0}, 

= {(a,^) e [a-,a-e] x K : (a,j/3) G 

OO 

J=2 

Here 91 is the set of all the purely imaginary zeros of all polynomials P that belong to the family 

([T§. 

We make the following assumptions. 

(RO) r is continuous in both variables for any r G r; 


(Rl) For any r G c. 


lim sup 

|| a :||->0 a 


lk(Q:,3^)11 
11*11 


= 0 ; 


(R2) For any a G [a_,a+], 0 ^ [ao(a), 5o(a)]; 

(R3) 0(a_) is ( 7 i-unstable (cf. Definition |2 .2 1 and ([ni)-®); 









(R4) i3(a+) is g 2 -unstable with qi ^ q 2 ] 

(R5') There exists a finite collection of disjoint sets Vk C [a_,a+] x K_|_ such that: 

(i) each Vk is homeomorphic to a closed disk; 

(ii) 

(iii) &(£/) Pi dVk = 0 for any fc; 

(R5") IK(^) n 6(^) = 0; 

(R5'") For any a € [a_, a+] and for any P G Q, P has at most one pair of purely imaginary roots. 

We are now in a position to formulate our main result on the Hopf bifurcation in interval 
systems. 


Theorem 3.7. Suppose that (R0)-(R4) hold and either (R5'), (R5'') or (R5'") is satisfied. 
Then, any selector (12) belonging to © has a branch of non-constant periodic solutions bifurcating 
from the trivial solution. 


Remark 3.8. One can easily see the parallelism between hypotheses of Theorem |3.7| and their 
non-interval counter-parts from Theorem |3.5[a,b,c). 


Remark 3.9. Conditions (R3) and (R4) can be verified using Kharitonov’s theorem (see Theorem 


2.3) and Lemma 2.4 To verify (R5'"), one can use Lemma 2.6 


4 Examples 


Below we present three examples illustrating Theorem 3.7 with one of the conditions (R5') - (R5'") 


in each of them. To simplify the exposition, we are dealing with higher order scalar equations rather 
than with equivalent first order systems. The class of nonlinearities r in each example is assumed 
to satisfy conditions (RO) and (Rl). 


The characteristic equation of the linearization of (16) at zero has the form 

Q(a)(A) = Jo(a) + Ji(a)A + J2{a)X^ + J3{a)\^ + AP 


Following (14), we compute 


Im(( 7 i(iJ(Q!), i/3)) = 4 — — (5 — 3a -I- -I- e)/3^ -|- /3^; 

Im(g 2 (Q(Q:)) i/d)) = 4 — 4a^ — (5 — 3a -I- a^ — e)/3^ -I- /3^; 

Re(/ii(Q(a), i/3)) = (4a — 3a^ -I- a^ — s)/3 — (a -I- a‘^)/3^] 

Re(/i2(0(a), i/3)) = (4a — 3a^ -I- a^ -I- e)/3 — (a -I- a^)/3^. 


(15) 


Example 4.1 (Theorem 3.7 with (R5')). Fix e = 0.28 and, for any real a, define four intervals 
as follows: 

Jo(a) = {4-4a^}; 

Ji (a) = [4a — 3a^ -|- a^ — e, 4a — 3a^ -I- a^ -I- e]; 

J 2 (a) = [5 — 3a -I- a^ — £, 5 — 3a -I- a^ -I- e]; 

•^ 3 ( 0 :) = {a -I- a^}. 

Consider the following forth order interval differential equation 

Jo{<y)y + Ji{a)y' + J 2 {a)y" + J 3 {a)y"' + y'"' = r(a, y, y', y", y'"). 


(16) 


(17) 
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Figure 2: (a) The dark grey domain that consists of two connected components is the set fH 

of purely imaginary characteristic roots i/3 of the interval polynomial ( |17[ ) for Example 4.1 The 
two solid black curves inside the two components of d\ show the set of purely imaginary roots 
for a representative polynomial P{a){ ) that belongs to the family 0 - This representative has 
two purely imaginary roots i/3i(a), i/ 32 (a) for some interval of a values [ 01 , 02 ] C (a_,o+) = 
(—0.4, 0.8). The light grey domains are the sets ©2 and 63 obtained from the dark grey domain 
91 by the transformations (o,/3) i-)- (o,/3/2) and (o,/3) !->■ {a, 13/3), respectively; the dashed curves 
inside ©2 are the images of the solid black curves in 91 under this transformation. The intersection 
of the solid curve and the dashed curve inside the smaller component of 91 corresponds to the 
2 : 1 resonance i/3i(o) = 2 // 32 (o). The dashed quadrangle Vi contains the set 91; its boundary 
does not intersect ©^ in accordance with (R5'). (b) The real parts 'ri(a), T 2 (q;) of the roots of 
the representative polynomial P{a,-) that belongs to the family 0 (schematic). The sliding 
intervals Ti{a) = 0, T 2 {a) = 0 correspond to the black curves (3i{a), / 32 (a) shown inside the dark 
grey domain 91 on panel (a). 


Take a_ = —0.4 and a+ = 0.8 and define 91 and © as in (14). Let us show that equation 


satisfies conditions of Theorem 3.7 with (R5'). Since, by construction, a € [—0.4, 0.8], (R2) is 

Observe that 


2.4 


satisfied (cf. the first formula in (15)). To show (R3) and (R4), we use Lemma^ _ 

i3(a)(A) from ( |T7| is obtained from the polynomial Pq( a)(A) = (A^+aA+l+a)-(A^+a^A+4—4a) by 
taking e-neighborhoods of some of its coefficients. By direct verification, Po(a+) is Hurwitz stable 
while Po(a_) is 2-unstable. To complete the verification of condition (R3) (resp., (R4)), it remains 
to observe that {(a_,/3) S [a_,a+] x M+j n 91 = 0 (resp., {(a+,/3) G [a_, a+] x K+j n 91 = 0). 
These last two relations as well as condition (R5') are illustrated by Figure]^ 


Example 4.2 (Theorem |3.7| with (R5")). Fix s = 0.7 and, for any real a, define four intervals 
follows: 

(/o(q^) = [81 T 27o; T — S', 81 T 27o; T 2o/^ T sjj 

(/i (q/) = [9o/ T 1 lo/^ “h — s, 9 gi T WoP" Ts]^ 

J 2 (g) = {18 + 3g + 

= [gTg — s,G-hG T s]. 


as 


As in Example 4.1 consider the interval differential equation (16) and the characteristic polynomial 
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(b) 


W2 


Wi 


Figure 3: (a) The nonintersecting sets 93 (dark grey) and © (light grey) for Example 4.2 with 
[Q!i,a 2 ] = [—0.5,1.2]. The black curve is the set of purely imaginary roots of a representative 
polynomial P(a)(-) that belongs to the family 0- At the corner point of this curve, has 

a purely imaginary root of multiplicity 2. The real parts of the roots of P(q;)(-) behave as shown 
on panel (b) of Figurej^ (b) Curves (wi(a±, ■),W 2 {a±, •)) (thick lines) for Example |4.3| The thin 
curve (wi(a, •), r(; 2 (a, •)) with a = 0.075 from the interior of the interval [a_,a+] intersects the 
negative cone {(■u;i,W 2 ) ■ wi < 0, W 2 < 0}. 


of its linearization ( [T7| ). In this case, 

Re(gi(Q(a),i/3)) = (81 + 27a + 2 ^^ _ e) _ (ig + 3a + a^)/?^ + 

Re{g 2 iQ{a), i/3)) = (81 + 27a + 2a^ + e) - (18 + 3a + a^)^^ ^ 4 . 

Im(/ii(iJ(a), i/3)) = (9a + lla^ + a^ — e)/3 — (a + a^ + e)/3^; 

Im(/i 2 (il(a), i/3)) = (9a + lla^ + a^ + e)/3 — (a + a^ — e)/3^. 

Consider the interval [a_, a+] = [—0.5,1.2]. For this interval, conditions (R2) - (R4) of Theorem 
|3.7| can be verified in the same way as in the previous example. In particular, one can use the 
representative polynomial Po(a)(A) = (9 + a + aA + A^)(9 + 2a + a^A + A^) when proving (R3) 
and (R4). Figure]^ shows that condition (R5") is also satisfied. 

Example 4.3 (Theorem |3 . 7| with (R5"')). Fix e = 1. For any real a, define five intervals 

Jo (a) = [36 — e, 36 + e]; 

Ji (a) = [36 + 36a — e, 36 + 36a + e]; 

J 2 {c() = [47 + 36a — e, 47 + 36a + ej; (18) 

J 3 (a) = [37 T 11a — £, 37 T 11a T S'jj 
J 4 (a) = [11 T a — £, 11 T a T S'] 

and consider the fifth order interval differential equation 

Ma)y + Ma)y' + J 2 (a)y" + Ma)y"' + Ma)y"" + y”'" = x{a,y,y',y\y"',y'"'). 

The corresponding characteristic polynomial equals 

Sj(a)(A) = Jo(a) + Ji(a)A + J 2 (a)A^ + J 3 (a)A^ + J 4 (a)A^ + A^. (10) 
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Let us take a_ = —0.1, a+ 
satisfied. By construction, e 
(R5'"), we apply Lemma |2 . 6 


= 0.09 and show that conditions of Theorem |3.7| with (R5"') are 
= 1, hence (R2) holds (cf. the first formula in (|l^). To show 
To this end, put Q(A) = A and i?(A) = —5. By direct calculation, 


r(Q, Q, R){uj) = (-144 - 180a ± 6 £)a; + (138 + 19a ± + (6 + a ± e)a;^ 


where for brevity we denote the interval [p, — £,/r + e] by /r ± e. Since for a G [a_,a+], T has 
at most one sign change, property (R5'") is satisfied. Finally, to show (R3) and (R4), we use 
the same argument as in the previous examples observing that 0 in (19) is obtained from the 
polynomial Po(a)(A) = (2 + A) • (3 + A) • (6 + A) • (1 + aA + A^) by taking £-neighborhoods of some 
of its coefficients. By direct verification, Po(a+)(-) is Hurwitz stable while Po(a-){-) is 2-unstable. 
To complete the verification of condition (R3) (resp. (R4)), it remains to observe that the c urve s 
shown in Figure]^ don’t intersect the negative cone {(a;, j/) G : x < 0, y < 0} (cf. Lemma 2.4). 


5 Proof of Theorem 3.2 


5.1 Necessary condition for the Hopf bifurcation 


Before proving Theorem |3.2[ let us show that the assumptions of Theorem 3.2 imply the classical 
necessary condition for the Hopf bifurcation. 


Proposition 5.1. Under the assumptions of Theorem 3.2. there exists (ao,/3o) G such that 

A(ao,/3o,0) = 0. 


Proof: Observe (cf. property (P4)(i)) that dV is homeomorphic to a 2-dimensional sphere. Take 

the standard orientation on and induce an orientation on Vo C dV. This orientation canonically 
induces orientations on V± and the orientation on dV. In particular, the local Brouwer degree 
for h\dv '■ dV —>■ C is correctly defined (provided that, say, the standard orientaion on C ~ is 
chosen). Since dV is compact and C is not compact, it follows that K\dv is not surjective and, 
therefore, 

deg(A,aP)=0 (20) 

(cf. [^, Chapter Vlll, Subsection 4.5). Combining ( [20| ) with condition (P4)(ii) and the excision 
property of the local Brouwer degree, one has (cf. [^, p. 277): 

deg(A,aP) = deg(A,P+) -hdeg(A,P_) -hdeg(A,Po) = 0. ( 21 ) 


By construction, the orientation on V+ (resp., V-) coincides with the orientation on Vo (resp., is 
opposite to it). Denote by i± the number of roots of A(a±, /3, r) in V± (counted according to their 
multiplicities). It is easy to see that i± = ±deg(A,P±). This observation together with formula 
( 21 ) implies 

deg(A,'Po) = t- - i+ 7^ 0 


(cf. condition (P4(iii)). On then other hand, combining condition (P5)(ii) with the Z 2 -equivariance 
of A (see condition (P4)) yields 


= ^deg(A,Po) = deg(A,|JX>fc). 


( 22 ) 


By the existence property of the Brouwer degree, the conclusion follows. 


□ 
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5.2 Normalization of the period 

We are looking for periodic solutions, with unknown period p, of the differential equation 

X = f{a,x). 


Following the standard scheme, let us introduce the unknown period p as an additional parameter. 
Define (3 = ^ and apply the change of variables 



to obtain the system 


u(0) = u{2tt). 


(23) 


We are now in a position to reformulate the original problem as an operator equation in the 
appropriate space of 27r-periodic functions and apply the equivariant degree method. 


5.3 S'^-representations 


We will use the first Sobolev space of functions on the unit circle equipped with the natural structure 
of S'^-representation induced by the shift in time. Let us recall some standard facts related to 5”^- 
representations. As is well-known (see, for example, 12 ), any real irreducible 
is of dimension 1 or 2 and can be described as follows. 


representation 
0 and define the 


S'^-action on C ~ by (e^^jz) >->■ e® • z, where 


Take an integer I > 
stands for complex multiplication, (denote 


this representation V;); also, denote by Vo the trivial one-dimensional 5'^-representation. 

Define V = K". Denote by W = the first Sobolev space of functions from to V. 

Observe that W admits the “Fourier decomposition” 


OO 

W = V©0Wi, 

Z=1 


(24) 


where the subspace of zero Fourier modes (i.e., constant functions) is identified with V, while the 
subspace of the Z-th Fourier modes Wi is identified with the complexification of V (denoted V^). 
In particular, any function u G Wi can be written in the form • {xi + iyi) for some xi,yi G V. 
There is a natural orthogonal S'^-representation on W given by 


(e^^,u)(t)<-^u(t + <p), e^^GS\uGW. 


(25) 


Formula (25) gives rise to the trivial action on V and the action (e*‘^,u)(t) i—>■ • u(t) on Wi. 


5.4 Reformulation in the functional space 

Take the first Sobolev space W and define the orthogonal projector K :W ^ L‘^{S^-,V) by 


1 

K{u) = — / u{t)dt, u G W. 
271 " Jo 


We can now rewrite (23) as the following operator equation in [a_, a+] x IR_|_ x W: 

/3,u) = u — 3F{a, f3,u) = 0, (a,/3) G [a_, a+] x K+, uGW, 


(26) 
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where 


T{a, j5,u) = {L + K) ^ ( — F{a,u) + Kii^ , 


L : W ^ L‘^{S^\ V) is given by L{x) = x and F : 


1 

X W 


(27) 


L‘^(S^;V) is defined by F(a,u)(t) = 
f(a,u{t)). Formula (251 gives rise to the S'^-action on [a_,a;+] x IR+ x W (we assume that 
acts trivially on [a_,a+] x M+). Moreover, it is easy to see that given by (26) and (27) is 


S'^-equivariant. 

5.5 Reducing the problem to computing S'^-degree 

In order to apply the equivariant degree method, we need to localize potential bifurcating branches 


in a cylindric box C [Q;_,a+] x IR_|_ x bF in such a way that the operator (26) is fl-admissible. 
To this end, consider the sets Vk provided by condition (P5) and put = dT>k. Since is 
compact, there are disjoint neighborhoods Nk of E*, C [a_,a+] x IR+ such that A(a,,0,O) ^ 0 for 
all (a,j3) G Nk- Set 


Vk=VkUNk-, ^k = dVk; V = \JVk; E = 91?; a(a,/3) = Id - Z?„.F(a,/3,0), (28) 

where D^F denotes the derivative of F with respect to u (cf. @). 

Lemma 5.2. There exists a disc Br{0) G W of radius r centered at the origin such that for all 
points (a,/3,u) G E x (Br(0) \ {0}), the following holds: 

(i) u — F(a, j5, u) 0; 

(ii) the fields ^(a,l3,-) and a{a, fi) are S^-equivariantly homotopic on Br{0). 


Proof: ^ 

(i) For a contradiction, suppose that for all p > 0, there exists {a,/3,u) G E x (Bp \ {0}) with 
u — F(a,l3,u) = 0. Since E is compact, without loss of generality, assume that there exists a 
sequence (aj, Pj,Uj) converging to (Q!*,,d*,0) such that (aj,j3j) G E, Uj — F(aj, f5j,Uj) = 0 and 
Uj 0. Then, 

PjFj) _ 

Ih.ll ■ WT"" ■ 

Observe that (aj,/dj) G E implies that /3j does not converge to 0. Combinig this with assumption 
(PI) and yields 


— Dj^F(aj, /3j, 0 ) ^ + 


r(aj,Pj,Uj) 


= 0 , 


(29) 


11“:; II ll“jll 11“:; II 

where r(aj, f5j,Uj) = F(aj, /3j,Uj) — DuF(aj, fij, 0)u. Also, 

U ' 'LL ' 

DuF(aj,l3j,0)^ = DuF(a,,/3,,0)^ + (7?„J-(a„/?„ 0) - D^F(a,, I3,,0))uj/\\u,\\. (30) 


Since DuF(a^, (3^,, 0) is compact, without loss of generality, we can assume that Z?„J^(a*, /?*, 0)uj/||Mj | 
converges to some u*. In addition, keeping in mind that Z?„J^(a*,/3*,0) depends continuously on 
a,/3, it follows from (30) that DuF(oij, l3j,Q)uj/\\uj\\ converges to u*. Combining this with ^ and 
(29) yields that converges to r;* 0. Hence (see ([^ once again). 


u* - Z?„J'(a*,/3*,0)u* = 0 
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meaning that Id — /3*, 0) is not invertible, which contradicts (P5)(iii). 

(ii) This part trivially follows from the compactness of S and condition (P5)(iii) combined 
with the standard linearization argument. □ 


Take given by (28) and 5^(0) provided by Lemma 


5.2 


Define 


Dfe = {(a,/3,M)G[a_,a+]xK+xIT:(a,/I)Gl?fc, weB,(0)}, D = |jDfc. (31) 

Clearly, Q, is S'^-invariant. By the existence of the invariant Urysohn function, one can take an 
invariant function c : D —>■ IR satisfying the properties 


c(a,/3,0)<0; i;{a,j3,u)>Q for ||m|| = r. 

Consider the map : D —>■ K 0 IT given by 

= (c(a,/3,M),g(a,/3,'u)). 


(32) 


By definition, any solution to the equation ^^(a,P,u) = 0 is also a solution to (23). In addition. 


is an S'^-equivariant D-admissible map for which 5'^-Deg(g,;,D) is correctly defined. 

Remark 5.3. As long as an invariant Urysohn function c satisfies properties (32), 5'^-Deg (g^, 0) 
is independent of the choice of c (homotopy property of the S'^-degree). 


The next statement provides a sufficient condition for the existence of a branch of periodic 


solutions bifurcating from the trivial solution (cf. Definition 2.1). We follow the scheme suggested 
in (see Theorem 9.18) with several modihcations making the argument more transparent. 

Proposition 5.4. Given s ystem assume conditions (PO) - (P5) are satisfied. Take fl defined 
by (31) and defined by (f5.5|)- Assume S'^-Deg (g,;, D) 7 ^ 0. Then, system has a branch of 


periodic solutions bifurcating from the trivial solution. 


As in [^, the following statement is the main topological ingredient in the proof of Proposition 5.4 
(cf. Theorem 3 in [^, p. 170). 

Proposition 5.5 (Kuratowski). Let X be a metric space, A,B C X two disjoint closed sets in 
X, and K a compact set in X such that KC]Af^tl)f^Kr\B. If the set K does not contain a 
connected component Kg such that Kg C A 7^ 0 7^ Kg C B, then there exist two disjoint open sets 
Vi, V 2 such that A C Ui , B C V 2 and AU B U K cViUV2- 

Proof of Proposition |5.4| Put 


ii: = g-i(o)nD. 

Consider the family of invariant functions Cq : D —>■ K given by 

i,q{a,f,u) = \\u\\-q, 0<q<r. 

Suppose for contradiction, there does no t exi st a compact connected set Kg C K with Aron<r(C^(0) 7 ^ 


KgGq^ ^(0). To apply Proposition 5.5 we need to show that AT n Cg ^(0) 7 ^ 0 7 ^ AT n ^(0). 


Not ice that for any q G (0,r), Cg satisf ies properties (32), so S'^-Deg (g^^, D 7 ^ 0 (cf Remark 
5.3 and the assumptions of Proposition 5.4 By the existence property of S'^-degree, for each 


q G (0,r) there exists {uq, j3q,Uq) G AT with ||ug|| = q. Since K is compact, it follows that there 
exist (ag, /Sg, Ug) G AT n C()'^(0) and (a^, fir, Ur) G AT n cuHO)- 
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Put 


A' = {{a,^,u) G n : ||■u||=0}, A” = {{a, f3,u) € ^ : ||u|| = r}. 

Then, by Proposition |5.5[ there exist open disjoint sets N' D A!, N" D A" with C iV' U N". 
Put 

Z’ = Z" = A” U (K n N') (33) 

and let us, first, show that Z' is S'^-invariant. Notice that A' is invariant. Suppose for contradiction 
that K C]N' \s not invariant. Then, there exist u € K n N' and 7 G S^ such that ( 7 , u) ^ K D N'. 
However, since K is invariant and K C N' U N”, it follows that ( 7 ,m) G K H N". We now have 
S^{u) C N'\J N" with S^{u) n ^ 0 ^ Si(u) n N" and N' n N" = 0, which contradicts the 
connectedness of S^{u). Thus, Z' is invariant as the union of invariant sets. Similarly, Z” is also 
invariant. 

Next, define an invariant Urysohn function /r : H > M with the following property: 


/i(a,/?,u) 


1, if (a,/3,u) G Z'] 
0, if (a, /3, u) G Z". 


(34) 


Take c : H 


defined by 


<,{a, /3, u) = ||u|| - ^(a, /3, u) ■ r. 

Clearly, c is invariant and satisfies properties 
tion 
with 


(35) 


(cf Remark 5.3 and the assumptions of Proposi- 
5.41 so S'^-Deg (l?^, fl) ^ 0. By the existence property of the S^-degree, there exists (a*, /3*, m*) 
— I3»,u^) = 0. Since iG C iV' U N” and N' n N” = 0, it follows that either 


(a*,/3*, u*)) G N' or (q!*,/3*,u*) G N". Assume (a*,/3*,ii*) G N'. Then (cf. (^l) and (^)), 
||m*|| = r, i.e. (a*,/3*,u*) G A" n = 0. Similarly, the assumption (a*,,0*,u*) G N" leads to a 
contradiction. 

5.6 Computation of S'^-Deg f2) via deformations 


Proposition 5.4 reduces the proof of Theorem|3.2|to the computation of S'^-Deg (1?^, fl) and showing 


that this degree is non-zero. Our goal now is to connect S'^-Deg (1?,^, O) to spectral properties of 
Dxf{ct,0) (cf. condition (PI)). This will be done in several steps. 


Step I: Reduction to a circle. Put U,;(q:,/ 3, u) = j3,u)^a{a, j3)u) (cf. (28l). Since <j(a,/?, 0) < 

’ii) that is iS^-equivariantly homotopic to on flk- 


0, it follows from Lemma 


5.2 


Take Di~ and Nk from (28) and assume, without loss of generality, that Nk is homeomorphic 


to [—1,1] X S^. Let (^fc, 7 fe) : Nk —>■ [—1,1] x 5'^ be a trivialization taking to {1} x S^. For any 
k, define three functions g*, : 2 ?^ —>• K, : flfc —>• K and ^ik : fife —)• M 0 W by 


9kia,l3) = 


0 , 


ia,l3)eVk\Nk-, 


(36) 

(37) 


^^kia,l3) + l, {a,(3)GNk; 

?k{a, P, u) = gkia, /?)(||'u|| - r) + ||u|| 0 r; 
di{a,l3,u) = (?fc(a,/3,u),a(a,/3)M). 

Obviously, the boundary dilk of the domain 0^ consists of three pieces: 

dUk = {||u|| = r} U {m = 0, ^kia,l3) = 1} U {u 7 ^ 0, ^kia, l3) = 1}. 

On the first piece, cj, and are both positive, while on the second piece they^are both negative. 
Also, on the third piece a{a,j3)u is non-zero. Hence, the vector fields and 5^,^ are not directed 
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oppositely on d^lk-, therefore they are equivariantly homotopic on 17^. Define x 5^(0). 

By (36) and (37), for all (a, j3) Nk, one has Qc > 0. Hence, by the excision and homotopy 

properties of the S'^-degree (see Appendix), S^-Deg (5^^, flk) = S'^-Deg 

Define rjk ■■ Nk ^ Efc by r]k{a,l3) = (^fc, 7fc)“H0> 7fc(a, ,0)) and -t K © IH by 


^k{a,/3,u) = {‘:k{a,X^,u),a{rjk{a,(3))u). 
Then, by the homotopy property of the 5'^-degree, 

S'^-Deg(S',;,Dfe) = S'^-Deg(^fc,Dfc). 


(38) 


(39) 


Observe that formulas (38), (39) reduce the computation of ^^-Deg (IJ,;, D^) to studying 5'^- 
equivariant homotopy properties of restrictions of Uk ■ Nk ^ GLf [W) to the zero section 
Sfe ~ {0} X where GL^ (W) stands for the group of S'^-equivariant linear completely con¬ 
tinuous vector fields in W. 


Step II: Computation of the degree. For any to € N, put bF™ = V (B Wi © ... © Wm (cf. (24)). 
Combining the compactness of the operator a with the suspension property of &e S'^-equivariant 
degree (see Appendix), one can find a sufficiently large to such that the field ^k is equivariantly 
homotopic to the compact field S'! : Dj, —>■ M © bF defined by 

^Ua,/3,u) = {gk{a,P,u),a{r]k{a,(3))u), 
where afe(a,/3) = afc(a,/3)|w™ + Id|(v^™)j_. Put 

aU<^,P) ■■=akia,/3)\v; ai{a, f3) :=akia, f3)\wi, I > 0. 

Fix some a between q;_ and a+. By condition (P3), the map aPj, : ^ GL{V) is homotopic 

to the constant map d : Nk ^ GL{V) given by d{a, f3) = —Dxf{d,0). Now, we are going to use 
formula (43) presented in Appendix. To this end, one needs to separate the “contribution” of the 
zero Fourier mode to the S'^-degree from other modes. Define 


nl = nlr\ © bF" 


ni = Bx nl, 


where B is the unit ball in V. Also, define —>• K © bF and —>• F © (K © bF) by 

771 

^k{a,(3,u) = {gk{a,P,u),^a[{r]k{a,P))u), 'Sl=dx ^k- 


1^1 


Combining the suspension property of the ^^-degree with the product formula (see [^, Theorem 
6 .8), one obtains 

^^-Deg Dfc) = S'^-Deg (S'fe, D^) = sign (det(a)) • S'^-Deg (S/c, ^D- 


Further, by applying formula (43), 

771 

^^-Degd?,;,^^,) = sign(det(a)) • ^ (deg(detc(ai),Dfc)(Z;)). 
Finally, applying the additivity property of 5^-Deg and the Brouwer degree, we get 

771 

^^-Deg D) = sign (det(a)) • EE (deg(detc(4),I?fe)(Zz)) 

k 1^1 

771 

= sign (det(a)) • E (deg(detc(ai),X>)(Zi)), 

1^1 


(40) 


17 









where ai(a,f5) = a{a, I5)\wi- 

Step III: Reduction to crossing numbers. Observe (see p. 266]) that 

ai{a,p) = -^[iipid- D^f{a,0)]. 

Put a'i{a,P) = il/3 ■ ai{a,(3). Since /3 > 0, the map a[ is homotopic to ai. Note (cf. §) that 
detc(aj(Q!,/?)) = A(a,/3,0). Finally (cf. condition (P4)(iii) and (22)), 

deg(detc(ai),:P) = ^ 0. 


Hence (cf. ^^-Deg (IJ.;, H) 7 ^ 0. The application of Proposition |5.4| completes the proof. 


6 Proof of Theorems 13.51 and 13.7 


6.1 Proof of Theorem 3.5K a,b,c) 


(a) Our goal is to construct a domain V satisfying (P4) in such a way that (P5') would imply 
(P5). To this end, take Vk provided by (P5') and a_, Q!+ provided by (P4'). Next, take a 
sufficiently large M > 0 to ensure that 


A-i(O) C{[a_,a+] xHm(O)} =: B, 


(41) 


where Bm{0) stands for the closed ball of radius M centered at the origin in the (/?, r)-plane. Also 
due to compactness, there exists a (5 > 0 such that 


A ^(0) n {0 < T < (5} C int((Jl?fe) X {0 < T < (5}. 


(42) 


Define 

V ={BC^{T> (5}) ulJPfc X {0 < r < (5}. 

Since B n Vk x {0 < r < (5} is homeomorphic to a disc, V satisfies (P4)(i). By the choice of 
M and 6 (see ( |^ and ([4^), V satisfies (P4)(ii). Also, (P4)' guarantees (P4)(iii). Finally, by 
construction, V and Vk satisfy (P5)(ii). 

(b) To prove Part (b), it suffices to deduce (P5') from (P5"). Notice that R{f) is the set of roots 
of polynomials with coefficients parameterized by a G [a_,a+]. Hence, the coefficients of these 
polynomials are uniformly bounded. Observe also that the leading coefficient of these polynomials 
is identically equal to 1, therefore R{f) is a compact set. 

For any £ > 0, there exists a sufficiently large m such that IJ^m ‘^*(/) O [«-,«+] x {/? < e}. 
Since [a_,a+] is compact and Dxf{-,0) is non-singular, it follows that R{f) is uniformly separated 
from [a_,a+] x {/3 < e} provided that e is small enough. On the other hand, the sets 
and i?(/) are compact and disjoint (see condition (P5")), so they can be uniformly separated. 
Hence there exists a neighborhood iVe(i?(/)) of R{f) in [a_, a+] x K_|_ such that iVe(i?(/))n5'(/) = 
0. Without loss of generality, one can assume that Ng{R{f)) is a finite union of discs, therefore, 
the complement to N^{R{f)) in [a_,a_|_] x M_|_ has finitely many bounded connected components, 
say, Set 

K 

P = 7V,(i?(/)) U U t/,. 

i-1 
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By construction, I? is a finite collection of disjoint sets homeomorphic to closed discs (denoted Vk) 
and dVi, C dNg{R{f)), thus Vk satisfies condition (P5'). Hence, the result follows from Theorem 

[Ha). 

(c) To prove part (c), it suffices to deduce (P5") from (P5'"). To this end, assume, by contradic¬ 
tion, that (P5") is not satisfied. Then, there exist a point {a, 13) S [a_,a+] x IR_|_ and an integer 
k >2 such that {a, (3), (a, /c/3) S R{f)- This contradicts (P5"'). 


6.2 Proof of Theorem |3.5K d) 

Let V be the set provided by condition (P4). Our first goal is to construct V D V such that (a) V 
satisfies (P4); and, (b) Vq is a disjoint union of finitely many sets homeomorphic to a closed disc 
(cf. §). To this end, without loss of generality (use a small perturbation of V if necessary), one 
can assume that Vq C [a_,a+] x K+ is a disjoint union where Ki is a (t'i -I- l)-connected 

compact domain. Using the same surgery argument as in the proof of Alexander’s tame sphere 
Theorem (see, for example, [^, Theorem 4.34), one can construct V satisfying (a) and (b). 

Our next goal is to construct a finite collection of discs Vk C [a_,Q;+] x K+ satisfying (P5). 
Take R and S given by (10). Using the same argument as in the proof of Theorem |3.5[ b) above, 
one can construct a sufficiently small neighborhood Ng{R Cl Vq) of the intersection RCiVq such 
that Ni;{R n Vq) n S = $ and Ng{R 0 Vq) C Vq (cf. condition (P5")). Take C = [a_, a+] x K+ \ 
Ns{R(^Vq). By the standard compactness argument, without loss of generality, assume that C 
splits into finitely many connected components C = Cq U Ci U ... U C^., where Cq stands for the 
(unique) unbounded component. Put V = U A^e(i?n Vq). Let us show that V is a. finite 

union of discs. By construction, V = u’^^iVk is a finite disjoint union of regular closed subsets 
(i.e., each Vk is a closure of its interior). To show that each Vk is contractible, take a closed 
curve 7 C I? and assume that it is not contractible to a point inside Vk- Then, there exists a set 
K C [«_, a+] X R^\'Dk bounded by 7 . However, this contradicts the construction of V. Therefore, 
Vk satisfies condition (P5)(i). Also, since dVk C 9A^E(i?n'Pg) for any k, it follows that Vk satisfies 
(P5)(iii). Finally, to show that Vk satisfies (P5)(ii), observe that RHV D RHVq. Since, V C Vq, 
one has RfiV = RHVq and (P5)(ii) follows. 


6.3 Proof of Theorem 13.71 


Clearly, if 0 ) satisfies (R0)-(R4), then any selector ( [T^ belonging to ( [Tl|) sa tisfies (P0)-(P4). 
Similarly, if satisfies (R5') (resp., (R5", (R5"')), then any selector ([l2[) belonging to ( [Tlj ) 
satisfies (P5') (resp., (P5", (P5"')). The result follows. 


7 Appendix: ^^-degree 

Let G be a compact Lie group acting on a metric space X (see, for example, [^). For any x G X, 
put G{x) = {gx G X : g G G} and call it the orbit of x. A set Z C X is called G-invariant 
(in short, invariant) if it contains all its orbits. Assume G acts on two metric spaces X and Y. 
A continuous map f : X ^ Y is called G-equivariant if f{gx) = gf{x) for ail x G X and g G G. 
In particular, if the action of G on U is trivial, then the equivariant map is called G-invariant. 
We refer to [^[^[^ (resp. [^ [I2|[^[^ ) for the equivariant topology (resp. representation theory) 
background frequently used in the present paper. 

Let V be an orthogonal S'^-representation. Suppose that an open bounded invariant set H C 
M 0 U is invariant with respect to the action, where we assume that acts trivially on M. We 
say that an equivariant map / : H —?> U is admissible if /“^{O} C dVi = 0. In this case, (/, H) is 
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called an admissible pair. Similarly, a continuous map : [0,1] x > 1/ is called an admissible 
(equivariant) homotopy if h{t,-) is admissible for any t G [0,1]. It is possible to axiomatically 
define a unique function S'^-Deg which assigns to each admissible pair a formal sum of finite cyclic 
groups with integer coefficients (cf. [^, pp. 109, 113 ). The following is a partial list of the axioms: 

(Al) (Existence) If S'^-Deg (/, fl) = nl^. ^ 0 for some k, then there exists an 

a; G O such that f{x) = 0 and Z 4 C Gx- 

(A2) (Homotopy) Suppose that /i : [0,1] x fl —>■ E is an admissible equivariant homotopy; then, 
S'^-Deg {h(t, •, •), n) = const. 


(A3)(Additivity) For two invariant open disjoint subsets C with / ^(0)nr2 C fli Ufl 2 , 

^i-Deg (/, n) = 51-Deg (/, fli) + 51-Deg (/, O2). 


(A4)(Normalization) Take Vi (cf. Subsection 5.3) and define the set flo and map b : IR 0 V 1 —>■ Vi 
by 


fin — 


|(t, 2 :)gK 0 Vi : |t| < 1 , 1/2 < ||z|| < 2|, b{t, z) = {1 - \\z\\ + it) ■ z. 


Then, S'^-Deg ( 6 , Oq) = 1 • (Zi). 


(A5) (Suspension) Suppose that A is an orthogonal S'^-representation and U is an open bounded 
invariant neighborhood of zero in A. Then, 


S'l-Deg (/ xld,nxU) = 5'1-Deg (/, fl). 


Using the equivariant version of the standard Leray-Schauder projection, one can define the S^- 
degree to 5'^-equivariant compact vector fields (see 5 26 for details). Combining the axioms of the 
S'^-degree with some standard homotopy theory techniques, one can reduce the computation of the 
S'^-degree of the maps naturally associated with the system undergoing the Hopf bifurcation to the 


computation of the Brouwer degree. To be more precise, let V be an orthogonal -representation 
with = {w G U : {j,v) = v V 7 G 5'^} = {0}. Take the isotypical decomposition 


u = Vfc, 0 Vfe, • • • 0 Vfc,, 


where each is modeled by the fc^-th irreducible representation. Define 

O = {(A,u) G C 0 U : ||u|| < 2 , 1/2 < |A| < 4}. 

Now, consider a map a : 5^ —> GL^^{V) and define aj : > GL^^ iVkj) by the formula aj{X) = 

a(A)n4^ (see, p. 284]). Let /q : O —>■ K 0 U be an ^^-equivariant map defined by 

/a(A,u) = ^|A|(||u|| - 1) 0 ||u|| 0 l,a . 

The following formula plays an important role in the proof of Theorem |3.2[ 

S 

S'^-Deg(/a,0) =^(deg(detcoaj,B))(ZfeJ, (43) 

1=1 

where B stands for the unit ball in C (cf. [^, Theorem 4.23). 
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